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Abstract 

ATHB17 (AT2G01430) is an Arabidopsis gene encoding a member of tlie a-subclass of tine liomeodomain leucine zipper class 
II (HD-Zip II) family of transcription factors. The ATHB17 monomer contains four domains common to all class II HD-Zip 
proteins: a putative repression domain adjacent to a homeodomain, leucine zipper, and carboxy terminal domain. However, 
it also possesses a unique N-terminus not present in other members of the family. In this study we demonstrate that the 
unique 73 amino acid N-terminus is involved in regulation of cellular localization of ATHBl 7. The ATHBl 7 protein is shown 
to function as a transcriptional repressor and an EAR-like motif is identified within the putative repression domain of 
ATHBl 7. Transformation of maize with an ATHBl 7 expression construct leads to the expression of ATHBl 7Al 1 3, a truncated 
protein lacking the first 113 amino acids which encodes a significant portion of the repression domain. Because 
ATHBl 7Al 13 lacks the repression domain, the protein cannot directly affect the transcription of its target genes. 
ATHBl 7Al 13 can homodimerize, form heterodimers with maize endogenous HD-Zip II proteins, and bind to target DNA 
sequences; thus, ATHB17A113 may interfere with HD-Zip II mediated transcriptional activity via a dominant negative 
mechanism. We provide evidence that maize HD-Zip II proteins function as transcriptional repressors and that ATHBl 7Al 13 
relieves this HD-Zip II mediated transcriptional repression activity. Expression of ATHBl 7A1 13 in maize leads to increased 
ear size at silking and, therefore, may enhance sink potential. We hypothesize that this phenotype could be a result of 
modulation of endogenous HD-Zip II pathways in maize. 
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Introduction 

HD-Zip proteins are a family of plant transcription factors that 
are found broadly across plant species [1] and that play an 
important role in regulating plant growth and development. In 
addition, many HD-Zip proteins are involved in mediating the 
response of plants to a range of environmental conditions. The 
family is characterized by the presence of a leucine zipper (LZ) 
domain adjacent to a homeodomain (HD) in the C-terminus [2,3] 



and comprises four distinct subfamilies designated as I, 11, III, and 
IV [1]. The sub-families, or classes, were defined based on the 
sequences of the conserved HD and LZ domains, gene structures, 
and the presence of additional class-specific motifs and domains. 
HD-Zip proteins bind DNA as homo- or heterodimers and their 
specificity of DNA recognition may depend on a precise spacing 
between the DNA-binding and dimerization regions [4]. Many 
HD-Zip proteins were shown to function as active repressors of 
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gene expression [4,5,6,7] and also to down regulate transcription 
of genes within the HD-Zip family [6,8] . 

Although the functions of many HD-Zip family members are 
poorly understood and their downstream targets are mainly 
unknown, recent studies have demonstrated a significant progres- 
sion toward the elucidation of potential fiinctions of HD-Zip 
proteins. The evidence indicates that these transcription factors 
play an important role in regulation of plant growth and 
development. HD-Zip IV proteins are expressed in the outer cell 
layers of various plant organs and are involved in the regulation of 
apical meristem and shoot organ development [9,10,11,12,13]. 
Members of the HD-Zip III family function in the maintenance 
and developmental regulation of apical mersteim and lateral organ 
polarity, most likely, through regulation of auxin transport 
[14,15,16,17,18,19]. Members of the HD-Zip I family have been 
associated with plant growth response to abiotic stresses 
[20,21,22]. The role of HD-Zip II proteins is still mainly 
undiscovered, as the null mutations in many HD-Zip II genes 
produce no detectable phenotypes [8,23,24]. Suggested functions 
for these proteins are based on studies of their expression patterns 
under dilferent conditions. ATHB2 was one of the first HD-Zip II 
genes described and shown to be regulated by the ratio of red to 
far red light (R/FR)[13]. ATHB4 is involved in die shade 
avoidance response [8]. Two other HD-Zip II proteins, HATl 
and HATS negatively regulate cell proliferation during leaf 
development under high R/FR light [25]. Additionally, HATl 
is also involved in reproductive development of Arabidopsis [26]. 
Additional functions of HATS and ATHB4 were elucidated using 
double mutations that led to simultaneous loss of both proteins, 
and it was demonstrated that they are involved in the regulation of 
cell proliferation during leaf development [27] and, potentially, 
regulation of shoot apical meristem activity and embryo develop- 
ment [28]. 

Systematic analysis by Zhao et al. [22] identified 18 HD-Zip lis 
out of 55 HD-Zip genes in the maize genome, but functional 
characterization of the maize HD-Zip lis has not been reported. 
While phylogenetic analysis suggests that HD-Zip II genes are 
closely related among monocot species, there are few ortholog 
pairs between monocot and dicot species [20,22]. Genes closely 
related to ATHB17 have been identified in rice and maize, 
OsHoxS and Zmhdz 1 8, respectively [22] . ATHB 1 7 is expressed in 
the root quiescent center of Arabidopsis [24,29], while OsHoxS is 
expressed in the panicle after flowering [20] . OsHoxS is reported 
to function as transcriptional repressors [30]. No physiological 
functions have been attributed to OsHoxS or Zmhdz 18 proteins 
and very littie is known about the potential functions of ATHB 1 7 
protein. No phenotypes were observed in a knoc-kout mutant of 
ATHB 17 by Hymus et al. [24]. In another study, seedUngs of an 
Arabidopsis knock-out mutant showed ABA-insensitive and 
drought-sensitive phenotypes, but only during a short post- 
germinating period, indicating ATHB 17 might function in a 
narrow window of plant development [SI]. Constitutive over- 
expression of ATHB 17 in Arabidopsis conferred several pheno- 
types including altered leaf development, changes in leaf 
morphology, plastid number and size, and increased photosyn- 
thetic capacity [24]. 

Here we demonstrate that ATHB 17 protein functions as a 
transcriptional repressor in Arabidopsis and corn protoplasts and 
contains a newly identified EAR-like motif responsible for 
repression activity. The protein is localized to both the cytoplasm 
and nuclei and its unique N-terminus could be responsible for 
regulation of protein localization. Similar to other members of 
HD-Zip II class, ATHB 17 forms homodimers and heterodimers 
within its own class and dimerization is required for high affinity 



DNA binding and repression activity. Transformation of the 
ATHB 17 gene into maize under the control of a monocot 
promoter results in the expression of a truncated protein that 
functions as a dominant-negative regulator that can modify 
activity of endogenous maize HD-Zip II transcription factors. 
Additionally, we demonstrate that maize HD-Zip II proteins 
function as repressors of gene expression and are expressed in a 
tissue- and development-specific manner. A greater number of 
HD-Zip II genes are transcribed in reproductive tissues, especially 
the c-ar inflorescence, than in any other tissues analyzed. 
Expression of ATHB17A1 IS protein in maize leads to changes 
in ear growth resulting in increased ear size at early reproductive 
stages and potentially increased sink size, which could be a result 
of ATHBl 7A1 13 modulating the activity of the endogenous maize 
HD-Zip II proteins. 

Results and Discussion 

ATHBl 7 is a Transcriptional Repressor Containing EAR- 
like Motif 

ATHB17 encodes a protein consisting of 275 amino acids 
(Figure 1) and contains five domains, four of which are topical for 
HD-Zip II proteins. Common to all HD-Zip II members is the 
presence of a homeodomain (HD) with an immediately adjacent 
leucine zipper (LZ), a hypothetical N-terminal repression domain, 
and a redox sensing motif (CPXCE), characteristic to all HD-Zip 
II proteins downstream of the leucine zipper motif [1,25]. In 
addition, ATHB 17 has a unique 73 amino acid N-terminal 
extension rich in tyrosine and cysteine that has not been identified 
in any other HD-Zip II proteins (Figure 1). 

To understand the potential functions of ATHB 17 protein, we 
attempted to characterize the unique N-terminal extension and 
determine if the protein has the ability activate or repress 
transcription. Based on bioinformatic analysis using the Trans- 
Mem-I- (GCG) program, the unique N-terminus of ATHB 17 may 
possess at least one, and possibly two, transmembrane domains at 
amino acids 4 to 23 and 35 to 57 or amino acids 4 to 23 and 40 to 
62, respectively. The presence of these putative transmembrane 
domains suggests that ATHB 17 could be tethered to a mc'mbrane 
[32] and translocated to the nucleus as a result of a specific 
regulatory mechanism [33]. Numerous examples exist of tran- 
scription factors that are sequestered in the cytoplasm, then 
tianslocated to the nucleus following a signaling event [32,34,35]. 
To determine if the N-terminus of ATHB 1 7 may play a role in 
mediating cytoplasmic or nuclear localization, C-terminal green 
fluorescent protein (GFP) fusions of ATHB 17 and a deletion 
variant lacking the unique N- terminus (ATHB17A73) were 
transformed into Arabidopsis protoplasts and the subcellular 
localization of each fusion protein was determined (Figures 2A 
and 2B). Deletion of amino acid residues 1 to 73, which contain 
the putative transmembrane domains, eliminated the cytoplasmic 
localization of ATHB 1 7 with only nuclear localization apparent. 
Two putative nuclear localization signals (NLS) were identified at 
amino acids 136 to 142 and amino acids 187 to 194 using the 
localization prediction tool pSort [36]. The putative NLSs were 
recognized as highly basic regions that is a common characteristic 
of NLS [3 7] . Arginines within each putative NLS were mutated to 
alanine and C-terminal GFP fusions of these variants were 
transformed into protoplasts. Only one variant (ATHB17A73 
R138A R142A) caused a loss of nuclear locaUzation (Figure 2C). 
Mutations of arginine to alanine at positions 187, 188, and 190 did 
not alter the nuclear localization of the A7S variant (Figure 2D). 
Therefore, the motif at amino acids 136 to 142 is the sole NLS 
necessary for targeting the ATHB 17 protein into the nucleus. 
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Figure 1. ATHB17 is a member of the o subclass within the HD-Zip II protein family. (A) represents the dendrogram and the domain 
architecture of the ATHB1 7 homologs. ATHB1 7 contains a typical homeodomain (HD; blue shading) and a leucine zipper motif (LZ; green shading) 
adjacent to the C-terminus of the HD. Red bars indicate conserved cysteines in the C-terminus. (B) shows the protein sequence of ATHB1 7. ATHB1 7 
contains a unique N-terminal extension (red shading) rich in cysteines and tyrosines. Additional structural feature identified for ATHB17 is a nuclear 
localization signal (red boxes). Downstream of the LZ motif is a putative redox sensing motif (CPXCE; red letters). 
doi:1 0.1 371/journal.pone.0094238.g001 



To determine if an active transcriptional repression domain is 
present in ATHB17, various deletions were constructed and tested 
in an Arabidopsis protoplast-based transcriptional activation/ 
repression assay. The pHATlv.GUS reporter construct was used 
because the basal level of expression allowed for observation of 
both activation and repression. Protoplast transformation with 
35S::chloramphenicol acetyltmnsferase {CAT) was used as a positive 
control because CAT is not expected to bind to DNA. 
Consequently, the amount of expression observed using CAT 
should reflect the basal level of reporter gene expression. 
Transformation of genes encoding proteins that bind to the 
promoter and cause activation or repression is expected to alter 
the basal reporter gene expression [38]. Results with specific 
eflFector constructs may vary between individual experiments from 
slight differences in transformation efficiency [38]. Transformation 
with ATHB2 was chosen as a positive control as it has been 
previously characterized as a transcriptional repressor [5,6,39]. 
Strong repression was observed when fuU-length ATHB 1 7 or fuU 
length ATHB2 was co-transformed with the /)/i4T7 ::&'[/>? reporter 
(Figure 3A). Expression of the pHATlv.GUS reporter was not 
repressed by ATHB 17 variants lacking amino acids 1-137 
(ATHB17 M-137) or 74-137 {ATHB17 h74-137) (Figure 3A), 
suggesting that all or part of a repression domain is present 
between amino acids 74 and 137. The increased expression 
observed in the presence of truncated ATHB 17 proteins would be 
expected if these truncated proteins compete for binding to the 
promoter with transcriptional repressors that are endogeiiously 
expressed in the protoplasts. 

The repression activity of several HD-Zip II proteins is 
associated with the presence of an ERF-associated amphiphUic 



repression (EAR) motif (LxLxL) [25], a well-described transcrip- 
tional repression motif [;40,41]. Two class II HD-Zips (ATHB2 
and ATHAT2) also have been functionally characterized as 
repressors using in vitro or in vivo assays [5,6,39]. However, an EAR 
motif was not identified in ATHB4, ATHATl 4, ATHB 17, or 
ATHB 18 proteins [25]. To characterize repression activity 
observed for ATHB 17 protein, we identified orthologous genes 
from several plant species and compiled sequences available in 
GenBank, including sequences for ATHB 17 and ATHB 18. A 
multiple sequence alignment of the N-terminus of these proteins is 
shown in Figure 3B. We were able to identify an EAR-like putative 
motif (LxLxI/M) at amino acids 84 to 86 related to the EAR 
transcriptional repression motif (LxLxL) [41]. A smaller putative 
EAR-like motif is located at amino acids 61 to 63. Mutating the 
conserved leucine residues L84 and L86 to alanine (L84A and 
L86A) caused a significant loss of repression (Figure 3C), while 
mutating leucine residues L61 and L63 to alanine (L61A and 
L63A) did not impact the level of repression (S. Tiwari, data not 
shown). This suggests that the EAR-like motif at amino acids 84 
and 86 is required to produce repression. 

ATHB17 Expression in Maize Results in a Truncated 
Protein that Lacks the Repression Domain 

The coding region of the full-length A THBl 7 gene under die 
control of a constitutive monocot promoter was transformed into 
maize. No signal corresponding to the fuU-length protein (expected 
molecular weight ~32 kDa) was observed in the transgenic maize 
plants on the Western blot shown in Figure 4. A band with an 
apparent molecular weight of approximately 20 kDa was clearly 
identified in the transgenic maize plants. The ~20 kDa signal 
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Figure 2. Representative fluorescence micrographs of Arabidopsis mesophyll protoplasts. Purified protoplasts were transformed with (Aj 
full length ATHB17:GFP or (B) a variant lacking the first 73 amino acids (ATHB17A73:GFP). The full length ATHB17:GFP fusion protein localizes to the 
nucleus and cytosol (A) in Arabidopsis protoplasts and ATHB17A73:GFP is localized exclusively in the nucleus when the unique N- terminus is 
removed (B). Two regions with a cluster of basic (K and R) residues were identified as putative nuclear localization signals. The arginine residues were 
mutated to alanines in these regions in the A73 variant to disrupt translocation. (C) R138A and R142A mutations in ATHB17 A73 disrupted nuclear 
retention while (D) R187A, R188A and R190A mutations did not alter the nuclear localization of the A73 variant. 
doi:10.1371/journal.pone.0094238.g002 



identified in tlie transgenic plants co-migrated witli the signal 
produced by the truncated ATHB17 protein produced in E. coli 
which was run on the same gel as a standard (Figure 4). As 
expected, no signal corresponding to the fiiU length or truncated 
ATHB17 proteins was detected in wild-type plants. A band with 
the apparent molecular weight of approximately 40 kDa was 
observed in both wild-type and transgenic events suggesting that 
the signal is due to the non-specific interaction between the 
ATHB17 antibody and endogenous maize proteins. Sequence 
analysis of the ATHB17 transcript produced in maize demonstrat- 
ed that, as a result of splicing, a truncated transcript was produced 
resulting in an ATHB 1 7 protein lacking the first 1 1 3 amino acids 
compared to the sequence of ATHB 17 found in Arabidopsis 
(hereafter referred to as ATHB17A113) (Figure SIA-SIC). 
Mapping of the maize ATHB 17 transcript sequence to the 
transformation vector revealed that the splicing donor site 
originates in the intron region of the construct, and that the 
splicing acceptor site is located within the coding region of the 
ATHB 17 gene, as shown in Figure SIC. Although in Arabidopsis 
no splice variants of the ATHB 17 gene were observed, the 



difiFerence in the recognition of the splice donor and acceptor sites 
between monocots and dicots is well described [42] . The truncated 
ATHB 1 7 protein identified by Western blot analysis corresponds 
to the size of the protein predicted to be translated from the 
transcript sequence (Figure SIB). 

As a result of the loss of the first 1 1 3 amino acid residues, the 
protein expressed in maize lacks its unique N-terminus containing 
putative transmembrane domains and a large part of the 
repression domain including the ERF-associated Amphiphilic 
Repression (EAR)-like motif, but has intact HD and LZ domains 
(Figure SID). Because ATHB17A113 retains dimerization and 
DNA binding domains but lacks the repression domain,we 
hypothesized that the truncated protein may affect the functions 
of the endogenous HD-Zip II proteins via a dominant-negative 
mechanism. To examine this hypothesis, we evaluated the 
properties of ATHB17A113 and the phenotypes associated with 
the expression of the truncated protein in maize. 
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Figure 3. Effect of various mutations on tKie transcriptional repression and identification of EAR-lil(e motif. (A) Arabidopsis mesophyll 
protoplasts were co-transformed with the pHAT1::GUS reporter gene and the following effector constructs: chloramphenicol acetyltransferase (CAT), 
ATHB17, ATHB17 with amino acids 74-137 deleted, ATHB17 with amino acids 1-137 deleted and ATHB2. Data are mean fluorescence readings 
measuring GUS-mediated substrate (4-methylumbelliferyl-beta-D-glucuronide) conversion. Error bars represent ± SD of three replicates. (B) Multiple 
alignment of the protein sequences similar to the amino terminus of ATHB17 from publically available HD-Zip II a subfamily proteins. Gaps are 
indicated by dashes. A putative EAR motif is shaded in red for the oc-subclass, while the EAR motif identified in [25] is shaded in green. Abbreviations 
for species are as follows: At Arabidopsis thaliana; Al Arabidopsis lyrata; Th Thellungiella halophila; Vv Vitis vinifera; Cp Carica papaya; Pt Populus 
trichocarpa; Rc Ricinus communis; Zm Zea mays; Sb Sorghum bicolor; Os Oryza sativa. (C) ATHB17, AJHBi? with the EAR motif mutated (ATHB1 7 L84A 
L86A), or CAT effector constructs were co-transformed into Arabidopsis mesophyll protoplasts with the pHATl::GUS reporter gene. Data are mean 
fluorescence readings measuring GUS-mediated substrate (4-methylumbelliferyl-beta-D-glucuronide) conversion. Error bars represent ± SD of three 
replicates. 

doi:1 0.1 371 /journal.pone.0094238.g003 
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Figure 4. Western blot analysis showing truncated ATHB17 
protein in transgenic maize events. Extracts from WT and 
transgenic maize events were separated on SDS-PAGE, transferred on 
to nitrocellulose membrane, and probed with ATHB1 7 antibody. An £ 
coli purified preparation of full-length and truncated protein was also 
analyzed for comparison with plant expressed protein as a control. 
Molecular weight markers are shown on the left. 
doi:l 0.1 371/journal.pone.0094238.g004 

ATHB17A113 Protein Localizes to the Nucleus in Maize, 
Binds to Consensus DNA Sequences, and Forms 
Homodimers and Heterodimers with Maize HD-Zip II 
Proteins 

To evaluate the localization of ATHB 1 7 A 1 1 3 in maize 
protoplasts, we transformed C-terminal GFP fusions of 
ATHB17A1 13 into maize protoplasts. A signal was observed only 
in nuclei, which suggests that, in planta, all expressed protein is 
transported into the nucleus (Figure 5). ATHB17A113 would not 
be exp(-( ted to be tethered to the membrane because of the 
absence of the putative transmembrane domains present at amino 
acids 4 to 23 and 35 to 57 or amino acids 4 to 23 and 40 to 62. To 
understand if the DNA-binding properties of ATHB17A1 13 have 
changed as a result of the truncation, we examined the DNA- 
binding properties of the ATHB17A113 in vitro using a surface 
plasmon resonance (SPR) assay. HD-Zip I and II proteins bind to 
similar cis elements (CAAT(N)ATTG) under in vitro conditions 
[4,30,43,44]. HD-Zip II members bind preferentially to the 
pseudo-palindromic sequence CAAT(C/G)ATTG [4], whereas 
members of the HD-Zip I subfamily bind preferentially to a 
pseudo-palindromic sequence that differs from that of HD-Zip II 
at the central nucleotide (CAAT(A/T)ATTG) [1]. As shown in 
Table 1, ATHB17A1 13 binds efficiently to both DNA targets with 
the equilibrium dissociation constant (Ko) of 37.7±14nM and 
20.4±3.3 nM for class I and class II DNA, respectively. The 
slightly higher affmit)' for class II DNA remained consistent 
throughout our kini;tic analyses, underscoring the class II DNA 
recognition preference by ATHB17A113. To determine if 
truncation of the N-terminus leads to changes in the specific 
sequence requirements for DNA binding, we examined the effects 
of mutating several nucleotides in the pseudo-palindromic 
sequence (Table SI). These data confirmed that ATHB17A113 
specifically recognizes the pseudo-paUndromic consensus sequence 
(CAAT(N)ATTG), consistent with previous characterization of 
oth(;r HD-Zip II proteins [4]. Additionally, we mutated amino 
acids within the HD (¥1820,185 N186) that should make specific 
base pair contacts with DNA [45]. As expected, mutating the 



residues corresponding to V182, 0,185, and N186 {ATHB17 
V182A Q185A m86A) completely abolished DNA binding 
(Table 1). 

The ability of ATHB17A113 to form homodimers was 
evaluated in maize protoplasts using a bead-based immunopre- 
cipitation assay with Luminex on-bead detection. A C-terminal 
MYC-HA fusion of ATHB 1 7 A 1 1 3 was co-transformed into maize 
leaf protoplasts with a C-terminal CFP fusion of ATHB17A1 13. 
The protein complexes were precipitated with anti-MYC or anti 
CFP magnetic beads and the expression of the protein of interest 
was confirmed with biotinylated antibody against HA or CFP 
(data not shown). Proteins in the dimer were precipitated with anti 
MYC magnetic beads and detected with l)iotinylated antibody 
against CFP. Transformation of either ATHB17A1 13::MYC-HA 
or ATHB17A1 13::CFP alone did not lead to a signal above 
background level. However, co-transformation of 
ATHB17A113::MYC-HA with ATHB17A1 13::CFP did lead to 
a signal (Figure 6A), demonstrating that ATHB17A113 can 
homodimerize. 

HD-Zip proteins have been shown to interact with members of 
their own family, but not with other families of HD-Zip proteins 
[4]. To identify maize proteins that potentially interact with 
ATHB 17, we performed a yeast two-hybrid screen of a maize 
cDNA library using ATHB17A1 13 as a bait. We identified several 
maize HD-Zip II proteins, but none of the ()tlii;r classes of HD-Zip 
proteins, as putative interactors with ATHB17A113 (data not 
shown). The results of the yeast two-hybrid screen indicate that 
ATHB17A113 may specifically interact with other HD-Zip II 
proteins when expressed in maize. To confirm this, we selected 13 
of the 18 HD-Zip II proteins as representative of the sub-classes in 
maize and tested them for interaction with ATHB 1 7 A 1 1 3 in 
maize protoplasts using the bead-based co-immunoprecipitation 
assay. 

The C-terminal MYC-HA fusion of ATHB17A113 was co- 
transformed into maize leaf protoplasts with a C-terminal CFP 
fusion of maize HD-Zip II coding sequences (Figure 6B). The 
proteins were exposed to antibody against MYC coupled to the 
magnetic beads and then precipitated. Antibody against CFP 
conjugated to biotin was used to determine if the complex 
included maize HD-Zip II proteins. As shown in Figure 3A, no 
signal above background level was identified when either 
ATHB17A113::MYG-HA or HD-Zip II::CFP was transformed 
alone. However, a clear signal was registered when 
ATHB17A1 13::MYG-HA was co-transformed with the constructs 
expressing HD-Zip II::CFP. No interaction was detected between 
HD-Zip I protein and ATHB17A113 (Figure 6B). These results 
confirm that when expressed in maize, ATHB17A113 can form 
heterodimers with maize HD-Zip II proteins and, therefore, has 
the potential to affect activities and pathways associated with 
maize HD-Zip II proteins. 

ATHB17A1 13 Lacks the Ability to Function as a 
Transcriptional Repressor 

Since ATHB17A113 lacks a large portion of the repression 
domain including the putative EAR-like motif, it should not be 
able to act as a transcriptional repressor. To determine whether 
ATHB17A113 protein has functional activit)' as a transcriptional 
repressor or activator in maize, we established a system that allows 
for detection of both repression and activation of reporter gene 
expression in maize mesophyll protoplasts. The transcriptional 
activity of ATHB17A113 was evaluated using two reporter 
constructs consisting of the GUS gene and either a Class I or a 
Class II binding sequence positioned between the 35S (—45) 
minimal promoter and the e35S double enhancer. A control 
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GFP cytosolic control 



CBF3::GFP nuclcarcontrol 



ATHB17A113::GFP 



Figure 5. ATHB17A1 13 is localized to nucleus. Fluorescence micrographs of maize leaf protoplasts transformed with (A) GFP alone, (B) AtCBF3 
fused to GFP as positive control for nuclear localization or (C) ATHB17A1 13 fused with GFP. Protoplasts from each construct were analyzed by laser 
scanning microscopy for GFP and chlorophyll fluorescence and acquired images were merged on bright-field image using the Image examiner (See 
IVlaterials and IVIethods for details). 
doi:1 0.1 371 /journal.pone.0094238.g005 



reporter construct contained neither a Class I nor a Class II 
binding sequence. Transcription was measured by co-transforming 
increasing amounts of ATHB17 expression plasmid with a 
constant amount of the reporter plasmid (Figure 7A). Repression 
of the class 11 promoter was observed with as little as 40 ng of 
ATHB17 expression plasmid and was dose-responsive up to 
630 ng of the ATHBl 7 expression plasmid. Repression of the class 
I promoter was observed with as little as 160 ng of ATHBl 7 
expression plasmid and increased at 630 ng. No repression of the 
control promoter lacking a binding site was observed. Thus, we 
confirmed that, in maize protoplasts, ATHBl 7 can repress 
transcription from promoters containing class II and, to a lesser 
degree, class I binding sites. Because no additional repression was 
observed above 630 ng of ATHBl 7 expression plasmid, we used 
smaller amounts in subsequent experiments to improve the assay 
sensitivity. 

In contrast to the repression of the GUS expression demon- 
strated by the fuU-length ATHBl 7 protein, no repression or 
activation of the ClassI::GUS or ClassILiGUS expression was 
observed when the ATHB17A113 plasmid was co-transformed 
with the reporter plasmids (Figure 7A). This result shows that the 
repression activity was lost as a result of splicing in maize. Based 
on these results, we conclude that ATHBl 7A1 13 lacks the ability 
to act as a transcriptional repressor. 



ATHB17A113 Acts as a Dominant Negative Regulator of 
Full-length ATHBl 7 Protein 

Although ATHBl 7A1 13 does not function as a repressor, the 
protein retains dimerization and DNA binding properties. 
Therefore, its likely action is to attenuate the activity of 
endogenous HD-Zip II proteins through a dominant-negative 
mechanism (Figure S2). The dominant-negative mechanism can 
occur through formation of non-functional homodimers or 
heterodimers with reduced DNA-binding activity or through 
competition for DNA-binding [46,47,48]. These mechanisms 
naturally occur and have been described in plants. For example, 
JAZ protein (JASMONATE ZIM-domain protein, a key repressor 
of Jasmonic acid signaling) has no DNA-binding domain, but can 
form heterodimers with the MYC2 transcription factor, a 
transcriptional activator. As a result, the complex between the 
JAZ and MYC2 proteins is transcriptionally inactive, and can 
repress the genes that are regxilated by MYC2 [49] . Formation of 
such competitive non-functional heterodimers with reduced DNA- 
binding affmity also has been described for Arabidopsis little 
zipper proteins, which can form heterodimers with HD-Zip III 
proteins. The heterodimers lack the ability to bind DNA and can 
reduce HD-Zip III protein activity [50]. Although formation of 
dimers with reduced DNA-binding activity is the most described 
mechanism thus far, accumulating evidence suggests that hetero- 
dimers containing DNA-binding domains but lacking activation 
domain(s) might reduce transcriptional activity as well as compete 
for DNA-binding with functional dimers [48] . 



Table 1. ATHB17An3 can bind both Class II and Class I DNA targets and ATHB17A113 containing V182A-Q185A-N186A mutation 
cannot bind Class II DNA target in in vitro assay (measured by Surface Plasmon Resonance (SPR). 





Protein 


Sequence 


ko„(1/Ms)x10' 


k„«(1/s)x10 " 


KofnlM) 


Relative Kp 


ATHBl 7A1 13 


CAGACA/^rOITTGCGGC (Class II) 


42.3 ±16.0 


880.6±416 


20.4±3.3 


1.0 


ATHBl 7A1 13 


CAGACA/irraTTGCGGC (Class 1) 


29.0±16.4 


1012.7±387 


37.7±14.0 


1.8 


ATHBl 7A1 13 


CAGCrCACTCTGAGGC (non-consensus) 






NB 




ATHBl 7A1 13 {V182A-Q185A- 
N186A) 


CAGACA/ATCATTGCGGC (Class II) 






NB 





Binding affinities and Kinetic constants of ATHBl 7A1 1 3 interacting with Class I and Class II type DNA, measured by Biacore 2000, globally fitted. SPR measurements witti 
Biacore 2000 were at 25°C in HBS-EP, 100 ug/ml BSA (10 mM HEPES pH 7.4, 150 mM NaCI, 3 mM EDTA, 0.005% Tween-20,100 ug/ml BSA). Equilibrium dissociation 
constant Kd = l<off/l<on. 
doi:1 0.1 371/journal.pone.0094238.t001 
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Figure 6. ATHB17A113 forms homodimers and makes heterodimers with maize HD-Zip II proteins. Maize leaf protoplasts were (A) 
mock- transformed or transformed with constructs expressing CFP, MYC::HA and CFP tagged ATHB17A113 alone or CFP and ATHB17A1 13::CFP in 
combination with ATHB17A1 13::MYC construct, (B) transformed (filled bars) with constructs expressing CFP or each of the CFP tagged HD-Zip lis 
alone or co-transformed (empty bars) with the construct expressing ATHB17A1 13::MYC-HA. Cellular extracts were co-immunoprecipitated using anti- 
MYC antibody and the immunoprecipitated complex was probed with anti-CFP antibody and visualized using phycoerythrin-conjugated reporter 
molecules. One HD-Zip I, Zmhdz3, was used to evaluate interaction with ATHB17A1 13. Nomenclature of HD-Zip II and HD-Zip I constructs is based on 
Zhao et al.[22]. 

doi:1 0.1 371 /journal.pone.0094238.g006 



To evEiluate the ability of ATHB17AI13 to act as a dominant 
negative regulator, we first examined the ability of ATHB17A113 
to relieve repression activity of the full-length ATHB 1 7 protein. 
We used the maize protoplast system described above to co- 
transform the reporter construct with 0.20 Hg of the fuU-length 
ATHB 17 plasmid and increasing amounts of ATHB17AI13 
plasmid (Figure 7B). Repression of the reporter gene expression 
caused by the full-length ATHB 17 was gradually relieved as 



increasing amounts of ATHB17A113 were added. This observa- 
tion suggests that ATHB17A113 has the ability to act as a de- 
repressor with a dominant-negative effect on expression of the 
reporter gene. This dominant-negative effect was only observed 
with the reporter containing a class II binding sequence, while no 
significant change in the level of expression was detected with the 
reporter containing class I binding sequence, or the reporter 
containing no class I or class II binding sequences. Since we 
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Figure 7. Full- length ATHB1 7 protein functions as transcriptional repressor and ATHB1 7A11 3 can relieve repression caused by full- 
length ATHB17 protein. Maize mesophyll protoplasts were transformed (A) with 4 |j.g cells of reporter (Class ll::GUS, Class l::GUS or No BS::GUS) 
and 0-5 |ig cells of effector (Full-length ATHB1 7) or 5 i^g of ATHB1 7A1 1 3 and Renilla luciferase (B) with 4 jig reporter (Class ll::GUS, Class l::GUS or No 
BS::GUS), 0-5 ^g ATHB1 7A11 3, and 0 (grey bars) or 0.2 yig (blue bars) of ATHB17full length. DNA amounts are per 320,000 cells. After 18 h, cells were 
assayed for GUS and luciferase expression. GUS values were divided by luciferase internal control values for each well and normalized to respective 
GFP samples. Bars are means and error bars represent 1 SD. 
doi:10.1371/journal.pone.0094238.g007 



observed lower affinity of binding to the class I sequence in the 
in vitro assay and an overall smaller effect on expression when the 
reporter construct contained a class I sequence, we hypothesize 
that ATHB17 may be most active against genes with class II 
binding sequences in the promoter regions. 

ATHB17A113 Acts as a Dominant Negative Regulator of 
Maize HD-Zip II Proteins 

Because ATHB17A1 13 can form heterodimers with maize HD- 
Zip II transcriptional factors and bind to the same DNA sequence, 
ATHB17A113 could be expected to act as a dominant negative 



regulator of endogenous maize HD-Zip II transcription factors. 
To evaluate this possibility, we cloned 1 3 members of the maize 
HD-Zip II family and tested their activity in leaf protoplasts as 
described above and in "Material and Methods". AU 13 HD-Zip 
II proteins evaluated in our study repressed GUS expression with 
both class I and class II binding sequences in a dosage dependent 
manner (Figure 8A). This observation suggests that maize HD-Zip 
lis act as a transcriptional repressor in vim. To determine whether 
ATHB17A113 has the ability to act as a dominant negative 
regulator with maize HD-Zip II proteins, the reporter constructs 
were co-transformed with 20 ng of each HD-Zip II construct and 
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an increasing amount of the ATHB17A113 expression construct. 
Dose-dependent relief of the HD-ZipII repression activity by 
increasing amounts of ATHB17A113 DNA was observed for all 
tested HD-Zip II proteins when the reporter contained a class II 
binding sequence in the promoter region (Figure 8B). In contrast, 
when the reporter construct contained a class I binding sequence 
or no binding sequences (control construct), no consistent dosage- 
dependent elfect of ATHB17A113 on the reporter expression was 
observed. We conclude that the ATHB17A1 13 protein can act as 
a dominant-negative regulator of endogenous maize HD-Zip II 
proteins that repress transcription of genes with promoters 
containing class II binding sites. This result further indicates that 
the mechanism of dominant-negative regulation exerted by 
ATHB17A113 should include a competition for class II DNA- 
binding sites. 

Expression of ATHB17A113 in Maize Leads to Increased 
Ear Weight at Early Reproductive Stages 

We investigated whether the expression of ATHB17A113 leads 
to any changes in maize plant development as HD-Zip II proteins 
are implicated in plant growth and reproductive development 
[26,51]. Field studies were conducted over two years in Illinois, 
USA. Standard agronomic practices (SAP) characteristic of corn 
production in the region were used to prepare and maintain plots. 
To determine whether ATHB17A113 affected plant growth and 
development in maize, several parameters were measured for two 
transgenic events in three hybrids. Because these are quantitative 
traits that may be subject to environmental variations, we analyzed 
all of our data across years and hybrids to maximize statistical 
power. As shown in Table 2, ther(; w(-r(- no significant changes in 
th(^ duration of time to tlu' rcprodu('ti\-c stage between transgenic 
events and their corresponding control, as evidenced by the 
number of days to anthesis and silking and calculated number of 
days between anthesis and silking (ASI). These results suggest that 
ATHB17A113 does not impact timing of maize development. 
There was, however, a significant increase in the ear weight at 
silking (the Rl development stage) in both transgenic events, as 
compared to the respective control plants (Table 3). There was no 
observed increase in either stover or total dry matter accumulation 
at Rl. The increase in ear weight without a concurrent increase in 
vegetative biomass also is evident in the increase in partitioning to 
the ear (Table 3). 

In maize, the sink size is determined by environmental and 
genetic factors around the Rl stage of development [52,53]. A 
tight genetic link has been suggested between sink size at flowering 
and final ^rain \'icld [54]. An increase in sink size and partitioning 
of assimilates to tlu' ("ar can k'ad to inc:reased yield if plant growth 
and development is supported by sulficient source during grain 
filling [55,56]. Environmental conditions during grain filling are 
critical for realization of yield potential [57,58,59,60], and 
contribute to the variability in measurements of yield. Hence, 
increased ear weight at Rl provides an opportunity for increased 
yield, l)ut the realization of this opportunity depends upon the 
environmental conditions during the grain filling period. 

Maize HD-Zip II Genes are Expressed in Ear Inflorescence 
and may be Involved in Growth Regulation 

Based on the molecular mechanism for ATHB17A113 
proposed above, it is likely that expression of ATHB17A113 
impacts maize ear growth by interfering with the functions of 
endogenous maize HD-Zip II proteins. Several HD-Zip II proteins 
have been implicated in regulation of growth and development of 
the meristematic tissue in Arabidopsis [61] including the recentiy 



identified JAIBA protein that was shown to regulate growth and 
development of the inflorescence meristem [26]. To im estigate 
whether the HD-Zip 11 proteins might be involved in regulation of 
ear growth in maize, we evaluated expression patterns of maize 
HD-Zip II genes in two different maize hybrids across several 
developmental stages in field-grown plants. The sampled tissues 
included leaf, leaf sheath, stem, ear inflorescence, tassels, ear, ear 
leaf sheath, cob, and kernels. 

The level of transcript in each tissue sampled was determined 
using the Quantigene 2.0 assay. The assay allowed for the 
simultaneous measurement of multiple RNA targets, and enabled 
analyses of the relative expression levels of each HD-Zip II gene. 
As shown in Figure 6A, maize HD-Zip II genes were found to be 
expressed in all tissue types and developmental stages similar to 
what was observed in rice [20] and Arabidopsis [13,22,25,26]. 
The number of HD-Zip II genes expressed in each tissue and the 
observed level of expression varied among HD-Zip II genes as 
shown in Figure 9A. Six genes (Zmhdzl9, 21, 25, 30, 33 and 35) 
were generally expressed in all tested tissues and developmental 
stages (Figure 9B), suggesting that they are expressed constitutively 
and, possibly, that they serve a more general role in the regulation 
of plant growth and development. Eight HD-Zip II genes 
(ZmhdzlB, 20, 23, 24, 26 29, 32 and 34) were predominantiy 
expressed in reproductive tissues (Figure 9C), suggesting that they 
might be actively involved in regulation of plant reproductive 
development. The maize ortholog of .\TIIB17, ZmhdzlS, was 
predominantly expressed in ear inflorescence and cob tissue, but 
not in kernels or the whole ear. This expression pattern is very 
similar to the expression of the rice ortholog, OsHox3, suggesting 
that the genes may have similar functions [20]. Based on their 
expression patterns, it is hypothesized that many maize HD-ZipII 
genes may play a role in developmental regulation during a 
specific period in plant development, as suggested for ATHB17 
[31]. 

Transcripts for a greater number of the HD-Zip II genes were 
detected in the ear inflorescence than in other maize tissues. 
Previously, we have shown in a protoplast system that all maize 
HD-Zip II proteins are repressors of transcription and that 
ATHB17A1 13 can modulate HD-Zip II functions and relieve the 
repression when co-expressed in the same tissue. Since 
ATHB17A113 is constitutively expressed in maize and its 
expression leads to an increase in ear size at silking, it is 
hypothesized that this phenotype could be a result of 
ATHB17A1 13 modulating the activity of maize HD-Zip II genes 
that are expressed in the ear inflorescence. Changes in the activity 
of maize HD-Zip lis could lead to changes in the expression levels 
of their targets. Although no direct targets other than HD-Zip lis 
were reported for over-exj)rcssion of full-length ATHB17 in 
Arabidopsis plants [31], ATHB17 may have small effects on the 
expression of its targets. Because ATHB17A113 only partiaUy 
inhibits the repression activity of endogenous maize HD-Zip II 
proteins, its effect on global transcription in maize would be 
expected to be very smaU. 

Expression of ATHB17A113 did not Lead to Major 
Changes in the Ear Inflorescence and Ear Transcriptome 

Because the ATHB17A1 13 expressing events have increased Rl 
ear biomass and because endogenous HD-zip II genes are aU 
expressed in developing ears, we hypothesized that ATHB17A1 13 
may influence the ear phenotype at Rl by modulating the activity 
of HD-zip II transcription factors. To evaluate this, we examined 
global transcriptome sequences in both ear inflorescence and early 
ear. We performed transcriptome sequencing using the lUumina 
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Figure 8. Maize HD-Zip II proteins function as transcriptional repressors and ATHB17A113 relieves the repression caused by 
endogenous maize HD-Zip II proteins. (A) Cells were co-transformed with 4 |ig/320,000 cells of each reporter (right side label) and an increasing 
amount of each HD-Zip II family member (left-to-right for each color: 0.008; 0.04; 0.2; 1; 5 |ig/320,000 cells). Error bars are 1 SD. (B) Cells were triple- 
transformed with 4 |ag/320,000 cells of each reporter (right side label), 20 ng of each HD-Zip family member and an increasing amount of 
ATHB1 7A1 1 3 (left-to-right for each color: 0.0, 0.3; 0.6; 1 .3; 2.5 ng/320,000 cells). For a given reporter, HD-Zip II, and concentration of ATHB1 7A1 1 3, the 
calculation of standard error relies on the number of biological replicates and the estimated error variance derived from a one-way ANOVA. 
doi:1 0.1 371 /journal.pone.0094238.g008 



HiSeq 2000 platform (Illumina, San Diego, CA) on two events 
in one germplasm. 

Each sample had between 1 2 million and 25 million reads with 
an average of about 18 milhon reads. Reads were mapped to the 
B73 public transcriptome (version 5a). On average, 29% of the 
reads mapped uniquely to a single transcript, 56% mapped to 
more than 1 but fewer than 1 1 transcripts, 1 % mapped to more 
than ten transcripts, and 14% did not map to any transcript. The 
reference transcriptome consisted of 136,770 transcripts. From the 
reference transcriptome, 76,612 transcripts were expressed (at least 



one count) in at least one sample, and 60,158 transcripts were not 
expressed in any sample. DESeq normalization scale factors 
ranged from 0.59 to 1.41. No samples were omitted due to quality 
concerns. As expected, many genes were differentially expressed 
between the V16 [62] ear inflorescence and Rl ear in the control 
plants (8961 transcripts; p<0.0001) and similar changes were 
observed for ATHB17A113 expressing transgenic events. How- 
ever, we observed very few transcripts that were differentially 
expressed between the ATHB17A1 13 expressing transgenic events 
and the control in either of the tissues. We identified 1 1 transcripts 
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Figure 9. Expression pattern of endogenous HD-Zip II transcripts in maize. (A) l-ieat map representation of endogenous HD-Zip II transcript 
expression across tissues and developmental stages in two maize hybrids. The red color gradient shows expression values represented on a Log2 
scale normalized median fluorescence intensity (MFI); a darker shade indicates a higher expression level. Grey color represents that expression was 
below background. For each of the listed tissue and development stages two hybrids (NH6214 and NN6306) are shown in the order that hybrid 
NN6306 is the first one. (B) Constitutive expression profiles of maize HD-Zip II genes. (C) Reproductive tissue specific profiles of maize HD-Zip II genes. 
Pink bar represents NN6306 and blue bar represents NH6214. Each panel shows the expression profile for the indicated HD-Zip II gene. The error bars 
are stderr. 

doi:1 0.1 371 /journal.pone.0094238.g009 



that were differentially expressed (raw_pS0.0001) in event 1, and 
24 transcripts that were differentially expressed in event 2 in ear 
inflorescence (Table 4). In ear tissue, 173 transcripts were 
differentially expressed in event 1 and five in event 2. To identify 
transcriptional changes that are most likely to be relevant to the 
phenotype, we looked for changes in transcript levels that were 
common to both events. Only seven transcripts showed differential 
expression (raw_ p^O.OOOl) in both events in the ear inflorescence 
(Table 5) and two transcripts showed differential expression in ear 
(Table 6). Functional assessment of transcript sequences by 
BLAST analysis indicates the differentially expressed genes in 
ear inflorescence have homology to plant Hsp22 heat shock 
protein, a protein containing RNA Recognition Motif (RRM), a 
BAG protein, and few hypothetical proteins. Small heat shock 
proteins can function as chaperones in vitro [63]. In maize, 
mitochondrial hsp22 has been suggested to protect mitochondria 
against heat stress [64] . RRM-containing proteins are involved in 
regulation of post-transcriptional gene expression [65]. BAG 
proteins are multifunctional and are hypothesized to be involved 
in cell regulation during growth, development, or stress response. 
Expression of BAG protein homologs was observed in a variety of 
tissues including developing flowers and ovaries in several plant 
species [66]. However, the biochemical pathways they may impact 
and specific functions of different family members are not yet well 
characterized in plants. Two transcripts were up-regulated in ear 
tissue of transgenic events (Table 6). One of the transcripts is 
annotated as glycerophosphodiester phosphodiesterase activity 
protein (GPPD). GPPD proteins are shown to play an important 
role in cell wall organization [67] . Cell wall organization is critical 
during cell expansion during growth. The second is annotated as a 
putative calmodulin-binding family protein which can be involved 
in transcriptional regulation [68] . It is not yet clear what role these 
proteins may play in ear growth and development in maize. 

Although HD-Zip II proteins have been shown to regulate the 
expression of each other [6,8], we did not observe any significant 
changes in transcript levels for any of the endogenous maize HD- 
Zip II genes in either tissue (raw_p<0.0001). Overall, the observed 
effects of ATHB17A113 on the maize ear inflorescence and ear 
transcriptome were very small. This is not unexpected, since the 
ATHB17A113 does not have direct repression activity and can 



only act through modulation of the activity of other HD-Zip lis 
and the overall effect could be very subtle. 

Conclusion 

ATHB17 protein functions as a transcriptional repressor and 
contains an EAR-like motif within the repression domain. Similar 
to other HD-Zip II proteins, ATHB17 requires a functional 
homeodomain and leucine zipper for transcriptional repression. 
The cellular localization of ATHB17 could be regulated through 
the N-terminus which is not present in other members of the 
family. Introduction of the ATHB17 coding sequence into maize 
leads to the expression of ATHB17A1 13, a truncated protein 
lacking a repression domain and transcriptional repression activity. 
ATHB17A113 forms homodimers which bind to the consensus 
DNA binding sequence, and heterodimerizes with endogenous 
maize HD-Zip II proteins. We show that the protein can disrupt 
HD-Zip II mediated transcriptional activity and is active in 
protoplasts as a dominant-negative regulator of the transcriptional 
repression activity of endogenous maize HD-Zip II proteins. 
Because ATHB17A1 13 lacks a functional repression domain, it is 
unlikely to have transcriptional activity by itself, and we 
hypothesize that ATHB17A113 may function by altering the 
activity of the endogenous HD-Zip II proteins. Modulation of the 
HD-Zip II regulated pathways may in turn lead to changes in the 
ear inflorescence growth and, ultimately, a larger ear at the early 
reproductive stage. The transcriptome analysis revealed only a few 
changes in the expression of native genes in the ear inflorescence 
and early ear, and no marked changes in HD-Zip II expression 
were observed. We therefore propose that subtie changes in 
transcription can affect agronomically important quantitative traits 
without causing undesirable phenotypic off-types that have been 
observed in past efforts to improve crop performance by 
overexpressing transcription factors [69]. 

Materials and Methods 

Ethics Statement 

Approval to conduct the field study in Unites States was 
obtained from United States Department of Agriculture. Field 
studies were performed on private land in Jerseyville (ILWl and 



Table 4. The number of differentially expressed genes in ear inflorescence and ear of ATHB17 events relative to control. 





Tissue 


No. of differentially expressed genes 










Event 1 


Event 2 


Overlap 


Ear Inflorescence 


Up 


4 


15 


1 




Down 


7 


9 


6 


Ear 


Up 


159 


3 


2 




Down 


14 


2 


0 



Table shows number of genes differentially expressed in individual event {raw p<0.0001) and common between events (overlap). 
doi:10.1371/journal.pone.0094238.t004 
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ILW2) and Wrights (ILWR), Illinois, United States of America 
(GPS co-ordinates for 2011-ILWl (39.0650, -90.3133), ILW2 
(39.0650, -90.3107), ILWR (39.3456, -90.2790) and for 2012- 
ILWl (39.0695, -90.3074), ILW2 (39.0695, -90.3074), ILWR 
(39.3381, —90.2856) for which Monsanto Company has lease 
agreement with land owner. Field studies did not involve any 
endangered or protected species. 

Comparative Sequence and Phylogenetic Analysis 

Protein sequences for ATHB17 orthologs and paralogs were 
identified by searching pubhc databases and were aligned using 
CLUSTALX [70]. Sequences belonging to the HD-Zip II a 
subclass (ATHB2, ATHATl and ATHAT2) were used as an 
outgroup. Transmembrane domains in ATHB17 were identified 
using TransMemH- (Wisconsin Package Version 10.3 software suite 
(Acct'lrys Inc., San Diego, CA, USA)) using the default settings. 

The protein sequence for ATHB 1 7 was analyzed for potential 
localization signals using pSort prediction software [36]. The 
software uses two rules to determine potential nuclear localization 
signals. The first is the identification of 4 residue pattern composed 
of basic amino acids (K or R), or composed of three basic amino 
acids (K or R) and H or P. The second is the identification of a 
bipartite signal composed of 2 basic residues, 10 residue spacer, 
and another basic region consisting of at least 3 basic residues out 
of 5 residues. The program identified a 7 residue pattern starting 
at amino acid 136. No bipartite signals were identified. A second 
potential nuclear targeting signal was identified by manually 
searching for highly basic regions. The region of amino acids 187— 
194 was identified composed of 5 basic residues out of 8. 

Transcriptional Repression Assay in Arabidopsis 
Protoplasts 

The 35S:: chloramphenicol cwetyltmnsferase {CAT) gene was described 
previously [38]. The cDNA sequence encoding fuU-length 
ATHB 17, ATHB 17 variants, or ATHB2 was cloned downstream 
of CaMF 555' promoter and a translational enhancer from the 5'- 
leader of Tobacco mosaic virus [71]. The promoter sequence for 
HATl (1821 bp) was amplified by PGR and was cloned upstream 
to the GUS [heta-glmuronidase) reporter gene. The pHATlv.GUS 
reporter gene was used in all transcription repression assays 
because the basal level of expression and magnitude of repression 
was better than that observed with pHAT2::GUS. The transcrip- 
tional terminator from the A. tumefaciens MOS gene was used in all 
reporter and effector plasmids [71]. The ATHB17 (AT2G01430) 
and ATHB4 (AT4G 16780) clones were PGR amplified from a 
diversified cDNA library created using mRNA from several 
different Arabidopsis tissues that was generated at Mendel 
Biotechnology, Inc. The ATHB17 PGR product amplified at 
Mendel contained an M227I substitution compared to the peptide 
sequence deposited at The Arabidopsis Information Resource 
(TAIR). This Mendel-derived sequence was used to construct all 
the variants described in this study. ^T//i?7 7 derivatives presented 
here were generated by standard techniques of molecular cloning 
[72] and the In-Fusion PGR cloning system (Clontech, Mountain 
View, GA). 

Protoplasts were prepared from whole leaves of three- to four- 
week olel Arabidopsis plants using an enzyme solution containing 
ceUulase (Research Products International Gorp, Mount Prospect, 
IL) and macerozyme (SERVA Electrophoresis GmbH, Islandia, 
NY) to remove cell walls [38]. For the transcriptional repression 
assay, 200 \)1 of protoplast cells were co-transformed with 10 |ig of 
the Gal4UAS[4X)::GUS reporter gene construct with either 10 |Xg 
of CAT, ATHB17, or ATHB17 variant constructs. Plasmid DNA 
used for co-transformation was prepared using EndoFree Plasmid 



Maxi kits (QIAGEN, Valencia, GA). DNA was introduced into 
protoplasts by incubating protoplasts with 40% PEG 4000 for 
30 min. After removal of PEG, protoplasts were kept in the dark at 

room temperature for 18 to 20 h. Protoplasts were then lysed in Ix 
cell culture lysis buffer (Promega, Madison, WI) and incubated 
with 1 mM 4-methylumbelliferyl-beta-D-glucuronide (Gold Bio- 
technology Inc., St. Louis, MO) solution at 37°G for 1 h. GUS 
activity was quantified by taking fluorescence measurements using 
a Synergy HT Microplate Reader (BioTek, Winooski, VT). Three 
biological replicates were performed and GUS measurements for 
each construct were averaged. 

Generation of Maize Transgenic Events and Expression 
Analysis 

A binary plasmid vector for transformation into maize was 
prepared containing the rice actin promoter, a leader and an 
intron driving the expression of fuU length ATHB17 coding 
sequence and 3'polyadenylation sequence. Primer sequences used 
to amplify ATHB17 coding sequence are as follows: ATHB 17 Fl- 
ATATATATggtaccATGATAAAAGTAGTATTTAGGTAG and 
ATHB17 Rl- ATATATATgggcccGGAGTGAAGGAT- 
GACGCTCTTG. Transgenic maize events expressing ATHB 17 
were developed through Agrobacterium-medmted transformation of 
immature maize embryos based on the method described by 
Sidorov and Duncan [73]. In brief, immature embryos were 
excised from a post-pollinated maize ear. After co-culturing the 
excised immature embryos with Agrobacterium carrying the ATHBl 7 
plasmid vector, the immature embryos were placed on selection 
medium containing glyphosate and carbenicillin disodium salt in 
order to inhibit the growth of untransformed plant cells and excess 
Agrobacterium. To confer tolerance to glyphosate the ATHB 17 
plasmid vector also contained the selectable marker GP4 EPSPS 
coding region regulated by rice actin promoter, leader, and intron, 
and a 3' polyadenylation sequence. Events containing the 
ATHBl 7 cassette were characterized by detailed molecular 
analyses to select events that contain only one copy of the 
ATHBl 7 cassette. 

To determine the expression of ATHBl 7 in transgenic events, 
proteins were extracted from the lyophilized leaf tissues from 
growth-chamber-grown plants using 1:50 extraction buffer con- 
taining 8 M urea. The samples were homogenized at 1200 rpm 
for 120 sec at room temperature (RT) using a megagrinder. They 
were then sonicated with a Model 2510 R-MTH Branson 
Ultrasonicator (Branson Ultrasonics, Danbury, GT) for 1 min at 
RT, centrifuged at 2800 rpm at 20°C for 20 min, and the 
supernatant was collected. Ten [i\ of sample were loaded into each 
lane of a 10% Bis-Tris LDS Midi gel (Invitrogen, Garlsbad, CA) at 
10 (il/lane as per MuPAGE Technical Guide. The proteins were 
transferred to a nitrocellulose membrane at 100 V constant 
voltage for 1 hour with IX transfer buffer containing 20% MeOH 
and 0. 1 % antioxidant. The membrane was blocked with 30 ml of 
WesternBreeze (Invitrogen) block for 2 hours at RT. The 
membrane was then washed twice with 40 ml ultrapure water, 
5 min each time, and probed with primary antibody diluted 
1:1000 in WesternBreeze block for 1 hour at RT with shaking. 
The membrane was washed thoroughly three times with 30 ml of 
antibody wash from the WesternBreeze kit for 5 min and probed 
with 30 ml of secondary antibody solution from the Wester- 
nBreeze kit for 1 hour at RT while shaking. The nitrocellulose 
membrane was then washed three times with 40 ml of antibody 
wash from the WesternBreeze kit for 5 min each and rinsed twice 
with 40 ml of ultrapure water for 5 min each. Detection was 
performed with the WesternBreeze GhemUuminescent Western 
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Blot Immunodetection Kit (Product#WB7106) as per the 
manufacturer's instruction. 

ATHB17 Transcript Sequencing 

Total RNA was extracted using a Trizol metliod. Briefly, Trizol 
reagent (Ambion by Life Technologies, Carlsbad, CA) was added 
to powdered leaf tissue and mixed by vortexing. Following an 
incubation at RT, chloroform (Fisher Scientific, Pittsburgh, PA) 
was added and mixed. Following another incubation at RT, the 
tube was centrifuged at 12,000xg at ~4°C for 10 min. The 
aqueous phase was transferred to a new tube, and the RNA was 
precipitated with a half-volume of isopropyl alcohol (Sigma, Saint 
Louis, MO). To pellet the RNA, the tube was centrifuged at 
12,000 xg at ~4°C for 20 min. The RNA pellet was washed using 
70% (v/v) ethanol and resuspended in DEPC-treated water 
(Ambion). P()ly(A) RNA was enriched using the Ambion Micro- 
Poly (A) Purist MAG Kit (Life Technologies, Carlsbad, CA) as per 
manufacturer's specifications. AH poly(A) enriched RNA was 
resuspended in water. Total RNA was quantified using a 
Nanodrop^'^ Spectrophotometer (Thermo Scientific, Wilmington, 
DE) following the manufacturer's instructions. cDNA was 
synthesized using Superscript III Reverse Transcriptase kit (Life 
Technologies, Carlsbad, CA) according to the manufacturer's 
specifications. 

A single PGR product (Product A) was generated that was 
analyzed to determine the nucleotide sequence of ATHBl 7 RNA. 
The PCR amplification of Product A was performed on event 2 
ATHB17 cDNA in a 50 ml reaction volume. Each reaction 
contained a final concentration of 0.2 |tM of each primer and IX 
Phusion HF MasterMix Buffer (NEB, Ipswich, MA). The following 
cycling conditions were followed: 35 cycles at 98°C for 30 sec; 
60°C for 30 sec; 72°C for 1 min with an initial denaturation step 
at 98°C for 30 sec and a final elongation cycle for 10 min at 72°C. 
Aliquots of each PCR product were separated on a 1% (w/v) 
agarose gel and visualized by ethidium br()mid(; staining to verify 
that the products were the expected size. Prior to sequencing, each 
verified PCR product was purified using an ExoSAP-IT PCR 
Product Cleanup Kit (ASymetrix, Santa Clara, CA) according to 
the manufacturer's instructions. 

Sequencing was performed using BigDye terminator chemistry 
(Applied Biosystems, Foster City, CA). A consensus sequence was 
generated by compiling sequences from multiple sequencing 
reactions performed on the PCR product. 

Localization of ATHBl 7 Full-length and Variants in 
Arabidopsis and Maize Protoplasts 

Arabidopsis mesophyll protoplasts were isolated [74] and 
transformed using polyethylene glycol (PEG4000). Approximately 
1.5x10'^ protoplasts were transformed with 10 mg of plasmid 
DNA per test construct for single protein expression. Maize leaf 
protoplasts were isolated from etiolated maize seedlings as 
described in Sheen and Bogorad [75] and transformed using 
PEG4000 [76]. Approximately 3x10^ protoplasts were trans- 
formed with 12 mg of plasmid DNA per test construct for single 
protein expression. Plasmid DNA was prepped using the Qiagen 
(Germantown, MB) Plasmid Mega kit. Protoplasts were incubated 
a total of 18 to 24 hours at 22°C and analyzed for GFP 
fluorescence using confocal microscopy. Confocal fluorescence 
microscopy analysis was performed using a Zeiss LSM5 1 0 META 
Laser Scanning Microscope (Carl Zeiss Microimaging, Inc., 
Thornwood, NY) equipped with a Krypton-Argon Ion (488 nm) 
laser, a green (543 nm) Helium-Neon laser, FITC, and Texas red 
filter sets. Image acquisition and analysis was performed using 
LSM 5 Image Examiner, version 4.2 (Carl Zeiss Microimaging, 



Inc.), using a 40x water 1.2 numerical aperture objective. 
Excitation wavelengths used were 488 nm and 543 nm, and 
emission filters were 500 to 530 nm and 630 to 700 nm to detect 
GFP and chlorophyll, respectively. 

DNA Binding 

Oligonucleotides, CAGACAATCATTGCGGC (Class II), CA- 
GACAATTATTGCGGC (Class I), and CAGCTCAGTCT- 

GACGGC (non consensus) were from Integrated DNA Technol- 
ogies (IDT, Coralville, lA). The pseudo-pahndromic sequences 
flanked by four nucleotides at the 5' end and 3 'end of the oligo are 
in bold. The Amine coupling Kit (BR- 1000-50) containing N- 
hydroxysuccinimide (NHS), l-ethyl-3-(3-dimethyl-amino-propyl)- 
carbodiimide Hydrochloride (EDC), Ethanolamine, pH 8.5, and 
Carboxylmethyl dextrane CMS sensor-chip, Ni[:kel-NTA (Ni- 
NTA) column, and HiLoad Superdex column were all from GE 
Healthcare Life Sciences (Piscataway, ]S[D, Streptavidin was from 
Thermo Fisher Scientific (Rockford, IL). Amine coupling of 
Streptavidin and capturing of biotinylated DNA were performed 
with 10 mM HEPES pH 7.4, 150 mM NaCl, 3 mM EDTA, and 
0.005% Tween-20 (HBS-EP). Biacore binding assays were in 
HBS-EP with 100 (Xg/ml BSA. Scrubber 2.0c was from Biologic 
Software Pty Limited (Campbell, Australia). Octet QK running 
buffer, IX PBS (10 mM Na2HP04, 1 mM KH2PO4, 137 mM 
NaCl, 2.7 mM KCl, pH 7.4) was from Roche Diagnostics 
(Indianapolis, IN). Tween-20 was from Sigma-Aldrich Corpora- 
tion (Saint-Louis, MO). Streptavidin sensor tips were from Forte 
Bio, Inc. (Menlo Park, CA) and microplates, black 96 well, were 
from E&K Scientific Products, Inc. (Santa Clara, CA). 

Tag free ATHBl 7A1 13 was expressed in E coli and purified by 
conventional chromatography. Surface plasmon resonance (SPR) 
experiments were performed using a Biacore 2000 instrument at 
25°C with filtered and degassed HBS-EP buffer. To generate a 
streptavidin sensor chip by amine coupling, CM5 sensor chip was 
conditioned with successive 20 sec injections of 0. 1 N HCl, 0. 1 N 
NaOH, and 0.5% SDS at 100 nl.min"'. The conditioned chip 
matrix was washed with HBS-EP buffer at a flow rate of 5 jtl 
min~', followed by activation of all four flow cells (Fc) with 35 |il 
of 1:1 NHS/EDC (75 mg mrVll.5 mg ml"'). Coupling was 
achieved with 35 |tl injection of 200 Hg/ml streptavidin in 10 mM 
sodium acetate (pH 6.5). Non-reacted carboxylmethyl groups were 
blocked with 35 nl of 1 M ethanolamine (pH 8.5). Biotinylated 
DNA (3 nM) in HBS-EP was captured at a flow rate of 5^1 min"\ 
The flow cell 1 (Fcl) was without captured DNA target and served 
as reference surface. 

Real time kinetics were performed at a flow rate of 50 |J,1 min~ 
in HBS-EP with 100 Hg/ml BSA as running buffer. Serial 
dilutions of five concentrations of ATHBl 7A1 13 (3.7 nM, 
11.1 nM, 33.3 nM, 100 nM and 300 nM) were injected. Two 
injections of buffer were performed before and aft(^r protein 
injections and used as blanks. Each 150 (tl protein/ huffc-r injc-ction 
was followed by a 400 sec dissociation period. The surface was 
regenerated for subsequent runs with a 10 nl injection of 0.01% 
SDS. All SPR data analyses were performed by global fitting to a 
1:1 Langmuir binding model using Scrubber 2.0c. 

All Octet-QK (Bio Layer Interferometry) experiments were 
performed in PBS (pH 7.4) supplemented with 0.005% Tween-20 
(PBS-T) and carried out at 30°C. Assay samples were in 200 ul per 
microplate well and agitated at 1 000 rpm. Streptavidin sensor tips 
were hydrated in reaction buffer before use. There are five steps in 
a typical Octet-QK assay. The initial three steps included a 
baseline step (1) to continue hydrating the sensor tips, followed by 
an immobilization step (2) to capture target DNA (6 nM) and an 
additional baseline step (3). Hydration, baselines and DNA capture 



PLOS ONE I www.plosone.org 



17 



April 2014 j Volume 9 | Issue 4 j e94238 



ATHBl 7 in Maize increases Ear Weight at Sill<ing 



steps were each for 10 min. In the last two steps which constitute 
the actual kinetic reaction, the association step (4) was initiated 
with samples of ATHB17A113 for 10 min and the dissociation 
step (5) was allowed for 15 min in PBS-T. Binding kinetics of 
target DNA variants were measured with 30 nM of 
ATHB17A1 13. Octet QK binding experiments were locally fitted 
to a 1:1 binding model using Octet 4.0 software package. 

Protein- protein Interaction 

Maize leaf protoplasts were isolated from 12- day-old plants as 
dcs(:rih(;d in Sheen and Bogorad [75] and transformed using PEG 
[76]. Approximately 3x10^ protoplasts were transformed with a 
total of 12 |ig of Qiagen prepped plasmid DNA 6 |a,g per 
expression construct or 6 |xg pGEM filler DNA. Protoplasts were 
incubated 18 to 24 hours at 22°C. Protoplasts were pelleted at 
150 g for 3 min and all but 300 |J,1 of incubation buffer was 
decanted. Pelleted cells were resuspended, 300 nl transferred to 
96-deep weU plate (#P-DW-500-C, VWR International, Radnor, 
PA) and centrifiaged at 150 xg for 3 min. The remaining 
incubation buffer was decanted, leaving 20 |a,l. The cells were 
lysed by adding 100 |tl Protein lysis buffer (50 mM Tris pH 7.8, 
150 mM NaCl, l%Triton X-100 (#9002-93-1, Sigma), Ix 
protease inhibitor Roche complete tablets (#11-836-153-001, 
Roche Diagnostics Corporation, Indianapolis, IN) and mixing 
three times every 15 min for 1 hour. The protein lysate was 
centrifuged at 3000 xg for 5 min and soluble fractions were 
retained for use in Luminex based co-immunoprecipitation (co-IP) 
assay (Luminex Corporation, Austin, Texas). 

Capture antibodies for CFP (D 153-3, MBL International 
Corporation, Woburn, MA), MYC (cat#A190-104A, Betiiyl 
Laboratories Inc., Montgomery, TX), and HA (cat#A190-107A, 
Bethyl Laboratories Inc.) were covalently coupled to carboxylated 
fluorescent microspheres (Luminex) according to the manufac- 
turer's protocol. Luminex two step Carbodiimide coupling 
protocol was modified as follows: 3.5 x 10^ beads were conjugated 
to 15 ng of antibody and resuspended in coupling buffer 
containing 100 mM MES pH 6. Biotinylated antibodies for 
CFP (#ab6658, Abeam), MYC (#A190-104B, Bethyl Laborato- 
ries Inc.) and HA (#A190-107B, Bethyl Laboratories Inc.) were 
used for detection of interacting prey proteins in the miniaturized 
sandwich immunoassay. NeutrAvidin R-phycoerythrin was used 
as reporter (A2660, Invitrogen). Protein expression for ATHB17- 
MYC::HA dual tag and CFP constructs were detected using the 
miniaturized sandwich immunoassay and co-IP method as 
described by Qi et al. [77]. Fifty |tl lysate was ahquoted to either 
3 or 4 different wells in a 96-well clear flat-bottom plate (#171- 
025001, Bio-Rad Hercules, CA). To each sample, 50 |tl of 
conjugated beads (5000 beads per well diluted in PBS-T 0.05% 
Tween and 0.5% BSA) was added to each appropriate well and 
incubated for 30 min with shaking. Forward and reverse Co-IPs 
were performed with conjugated beads mentioned above. Beads 
were captured for 2 min and washed 3 times with PBS-T. 
Biotinylated antibodies at (1:1000) were added to appropriate 
weUs. The plate was incubated on a shaker for 30 min followed 
by bead capture and wash. A 100 nl of reporter NeutrAvidin R- 
phycoerythrin at (1:1000) was added to each well and incubated 
30 min on a shaker with subserjuent washing. Each well was 
resuspended in 100|j1 of PBS-T buffer before plates were 
analyzed with FLEX MAP 3D system (Luminex Corporation). 
Approximately 100 beads were measured per sample to 
determine the median fluorescence intensity (MFI). 

Samples of mock transformed protoplasts, bait alone, and prey 
alone were used to determine the background signal of the 
Luminex assay. The mean and standard deviation of the MFI of 



these negative controls were calculated based on biological 
duplicates or triplicates. The mean+3 sigma was used as cut-off 
to distinguish background signal from positive signal. Signal above 
the +3 sigma background was designated as positive signal. 
Comparisons were performed using ANOVA t-test (JMP 8, SAS 
Institute, Gary, NC). 

Yeast-two-hybrid Screen 

ATHBl 7 full-length and multiple fragments were used as bait 
in a yeast-two-hybrid screen performed using ULTImate Y2H^'^ 
by Hybrigenics (Paris, France). cDNA encoding full-length 
ATHBl 7 was cloned into pB27 and pB29 (LexA- C and N 
terminal fusions), EAR motif region 1—91 into pB29 (LexA N- 
terminal fusion), Homeodomain region 128-234 into pB27 (LexA 
C -terminal fusion). Leucine zipper and Homeo-domain region 
114-275 into pB27 (LexA C-terminal fusion). These bait 
constructs were used to screen a random-primed cDNA library 
(MEGA_RP1) prepared from maize RNAs from callus tissue, 
etiolated seedlings, V3 seedlings, ear inflorescence, developing 
kernels, and ear leaf. 1 1 1-139 million clones (1 1-14 fold coverage 
of the library) were screened. cDNA fragments corresponding to 
positive 330 "prey" clones were amplified by PGR and sequenced 
at their 5' and 3' junctions. The resulting sequences were 
searched against Monsanto's proprietary database and assigned a 
quality score (A through F) which was indicative of the confidence 
of interaction. 

Transcriptional Repression Assay in Maize Protoplasts 

Maize protoplasts were transformed as described above. 53,000 
ceUs were transformed for each replicate. Each treatment was 
tested in four technical replicates per transformation. Biological 
replicates consisted of identical treatments tested in protoplasts 
isolated on different days. Cells were transformed with up to 6 
constructs: varying amounts of 1-2 effector DNA constructs 
(amounts are indicated in each figure), a constant amount of one 
GUS reporter construct, varying amounts of GFP, and constant 
amounts of two internal controls. GFP fluorescence was not used 
as a reporter in these assays; GFP was included to balance both the 
total amount of DNA and the amount of 35S-equivalent 
expression in each [:ell. Hence, when titrating effects of the 
various effectors, differences in amount of DNA transformed per 
well were normalized by co-transforming with GFP so that each 
well contained no more than 10 ng/320,000 cells of DNA. After 
addition of DNA, the total volume in each well was normalized by 
the addition of ddH20. Internal control constructs expressing 
firefly and Renilla luciferase, were used to normalize for normal 
well-to-well variations in transformation efficiency. Following 
background subtraction using "cells only" samples, data was 
normalized by dividing GUS counts per second per weU by 
relative luciferase units per weU. 

Physiological Phenotypes at R1 

Field experiments were planted in lUinois, USA, in 2011 and 
2012. Two different ATHBl 7A1 13 expressing events were 
compared to control plants in the context of three different 
Monsanto elite hybrids (NH6214, EXP257, NN6306). When 
insect protection and gK^)liosatc tolerance' traits were included, 
both test and control entries contained the same traits. NH6214 
and EXP257 were tested in 201 1 and NH62 14 and NN6306 were 
tested in 2012. The experiment was a randomized complete block 
design (RCBD) with 10 replicates and two testers blocked 
separately within the trial in 2011. The experiment was a group 
unbalanced block design (GUBD) with 18 replicates and two 
testers blocked randomly within the trial in 2012. Experimental 
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units consisted of three 2-rovv plots (21.5 ft long x 15 ft wide) 
blocked together for a total of 6 rows for each entry in the study. 
Standard agronomic practices (SAP), which are characteristic of 
corn production practices for the region, were used to prepare and 
maintain each field .site. Maintenance pesticides were apphed as 
needed and all maintenance operations were performed uniformly 
over the entire production area at a given field site. Plant 
development was followed on 10 tagged plants per plot. Days to 
reproductive developmental stages were assessed from planting 
until 50% of the plants reached anthesis and silking. Anthesis- 
silking interval was calculated by subtracting the number of days to 
50% anthesis from the number of days to 50% silking. 

Biomas.s samples were collected when plants reached the Rl 
growth stage. In 201 1, ten plants from a given area were harvested 
for biomass and the area measured; while in 2012, all plants from 
a 1-m row were sampled for biomass and the plant number was 
counted. Plants were harvested by cutting the stalk at soil level and 
separating each plant into its component parts, including leaf 
blades, stalks with leaf sheaths, and ear shoots (with husk and 
shank). AH plant parts were dried at approximately 70°C until a 
constant weight was achieved. Biomass of each plant component 
was recorded. Stover biomass was calculated by summing biomass 
of leaves and stalks. Total biomass was calculated by summing 
biomass of stover and ear shoots. Data from 201 1 was presented as 
g/m^, while data from 2012 was converted via covariate analysis 
into plant number adjusted g/m^. The ratio of ear weight to total 
weight was calculated at Rl to determine relative early 
partitioning of assimilates to the developing ear. 

Data were analyzed with the MIXED procedure of SAS 9.2 
software (SAS Institut(;, Gary NC). Outliers were identified using 
the deleted studentized residual (DSR) method, and outUer data 
points were discarded before further statistical analysis. Each data 
set consisted of events 1 and 2 and the corresponding control. 
Entries were analyzed across environments by using a mixed 
effects model. For model fitting, hybrid entry was considered a 
fixed effect, while environment, rep within environment, the 
environment x entry interaction and the residual (experimental 
error) were considered random elFects. Environment was defined 
as the combination of set within year and field. The statistical 
model used was: 

Yijr=U+ Li + S(L)„. + Ej + LEij + e,>, 

where Yyr = Observation from the j-th entry within the r-th rep at 
the i-th environment, U = overall mean effect, Ej = Effect of the j- 
th entry, L; = random effect of the i-th environment, B(L)ri = ran- 
dom effect of the r-th replicate within the i-th environment, LE 
ij = random interaction effect between the i-th environment and 
the j-th entry, and e = experimental error (residual) correspond- 
ing to the observation Yij^. 

Point estimates of the mean of event and control, difference 
between event and control means (delta), delta as a pcrcc-ntage of 
control mean (percent delta), and p-value for the test of delta were 
provided for comparing event vs. control across 2011-12 and 
across hybrids, and by year and by hybrid. 

HD-Zip II Expression Analysis by Quantigene 2.0 

Tissue samples were collected from 2012 field study (hybrids 
NH6214 and NN6306). RNA was extracted using the EZNA 
RNA Purification Kit (Omega BioTek, Norcross, GA). The 
transcript expression of 18 HD-Zip lis was analyzed using 
Quantigene 2.0. Raw mean fluorescence intensity (MFI) were 
obtained using the Quantigene 2.0. QC step included assigning 
0.5 LOD value for MFI that are below LOD only if more than 



50% of the data available. LOD was calculated using the formula 
LOD = AVG MFI of assay background control wells-l-3 SD of 
assay background signals. Each raw MFI was then subtracted with 
the averaged no template background MFI. 

This background subtracted MFI was normalized prior to log 
transformation. For normalization, the mean of housekeeping 
gene was computed across all samples (Sentries xSreplicates). Then 
a correction/normalization factor was obtained by dividing mean 
of housekeeping gene by each individual MFI of housekeeping 
gene (Asparaginase). This correction/normalization factor was 
multiplied with each individual background subtracted MFI for 
each "trait" (datapoint within that tissue type) and then Log2 
transformed to get the final normalized values for each trait. 

RNA Extraction, RNA-seq Sample Preparation and 
Sequencing 

Tissue samples were collected from the 2012 field study 
(hybrid NN6306). Total RNA was isolated from three biological 
repUcates for each tissue using TRIzol (Invitrogen) according to 
the manufacturer's protocol. RNA quantity was determined with 
a Nanodrop 8000 spectrophotometer and integrity assessed by 
Bioanalyzer assay with RIN greater than 6. Two |J,g of total 
RNA was used for sequencing library preparation with lUumina 
TruSeq RNA Sample Prep Kit V2 following manufacturer's 
protocol. qPCR (SYBR PGR Master Mix, Applied Biosystems) 
has been utilized to quantify sequencing libraries. Sequencing 
was performed with HiSeq2000 sequencing using 50 base pair 
reads. 

Sequence Alignment, Data Normalization, QC Analysis 

and Statistical Analysis 

The Bowtie apphcation was used for mapping the reads to the 
public B73 maize transcriptome reference. We allowed up to 2 
high quality mismatches for every 28 seed sequence, where "high 
quality" means an average Phred score greater than 20. Individual 
reads were mapped to the transcriptome. Reads may map 
uniquely, to multiple transcripts, or to no transcripts. If a read 
mapped to more than 10 transcripts, it was ignored. If it mapped 
to 10 or less, it was assigned equally to all of them. The DESeq 
algorithm (Simon Anders, EMBL) was then used to normalize the 
raw counts across all the samples in an experiment. Principal 
component analysis, hierarchical cluster analysis, and correlation 
analysis were used for QC measures and identify poor quality 
samples. Proc GLIMMIX (in SAS) was used to identify 
differentially expressed genes. 

The data discussed in this publication have been deposited in 
NGBI's Gene Expression Omnibus [78] and are accessible 
through GEO Series accession number GSE51742 (http://www. 
ncbi.nlm.nih.gov/ geo/ query/ acc.cgi?acc = GSE5 1 742). 

Supporting Information 

Figure SI (A) ATHB17 Transcript sequence from transgenic 
event. (B) Predicted Amino Acid Sequence of the Protein 
Produced by ATHB17 RNA in transgenic event. The generated 
consensus sequence of ATHB 1 7 transcript from the transgenic 

event \vas usc-d in a BLASTX 2.2.23 search of the GenBank_- 
Protein_Preferred database. The predicted amino acid sequence of 
ATHB 17 protein starting with first methionine is shown. 
Compared to the Arabidopsis ATHB 17 protein sequence 
(GenBank ID 179876107) , the ATHB 17 predicted protein 
sequence is truncated in transgenic event; lacking the amino- 
terminal 113 amino acids and is designated as ATHB17A1 13. (C) 
ATHB 17 transcription and translation products in transgenic 
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maize events. ATHBl 7 mRNA sequence identifies truncation of 
the 5' region of the ATHBl 7 coding sequence. This transcript 
sequence is predicted to produce a truncated protein lacking the 
first 1 13 amino adds. (D) Domain structure of ATHBl 7. ATHBl 7 
coding sequence contains canonical Homeodomain and Leucine 
Zipper (LZ) domains of HD-Zip family. Homeodomain is required 
for DNA binding. LZ domain is responsible for homodimerization 
and hetero-dimerrization with other HD-Zip II proteins. A 
repression domain is present upstream to the homeodomain. 
ATHBl 7 protein expressed in transgenic maize events lacks the 
first 1 1 3 amino acids resulting in deletion of part of repression 
domain containing EAR-like motif 
(TIF) 

Figure S2 ATHB17A113 will function via dominant 
negative mechanism. ATHBl 7 is expressed as truncated 
protein in maize lacking part of the repression domain. Maize 
endogenous HD-Zip lis function as transcriptional repressors. 
ATHB 1 7A 1 1 3 can interact with endogenous HD-Zip lis and 
sequester endogenous proteins from binding to their targets 
resulting in relief of repression caused by maize HD-Zip lis. In 
addition, heterodimer of ATHBl 7A1 13 and endogenous HD-Zip 
lis or ATHBl 7 Al 13 homodimers can compete for DNA binding 
resulting in altered target expression due to inability to cause active 
repression. 
(TIF) 

Table SI Evaluating ATHB17A113 binding affinities 
with various mutations in targets DNA measured by 
Octect^^K; locally fitted. 

(DOCX) 

Table S2 Phenology o{ ATHBl 7 events and control. The 

number of days to 5Q% silking and anthesis were m(;asured and 
the number of days between anthesis and silking ^^•as calculated 
(ASI) each year for physiological studies conducted under standard 
agronomic practices conditions. Differences in phenology between 
ATHB 17 events and control were determined using a fixed effects 
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